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WAVE REFRACTION DIAGRAMS FOR THE BALTIMORE CANYON REGION 
OF THE MID-ATLANTIC CONTINENTAL SHELF COMPUTED 
BY USING THREE BOTTOM TOPOGRAPHY 
APPROXIMATION TECHNIQUES 

Lament R. Poole 
Langley Research Center 

SUMMARY 

■ 

■ r ■ 

Wave rje/raction diagrams were computed by applying the 
Langley Research Center and Virginia Institute of Marine Science 
wave refraction computer model to the Baltimore Canyon region of 
the mid-Atlantic continental shelf. Diagrams were presented for a 
wide range of normally expected wave periods and directions. For 
each combination of period and direction, diagrams were computed 
by using each of three commonly used bottom topography approxima- 
tion techniques (quadratic least squares, cubic least squares, and 
constrained bicubic interpolation), and all these diagrams are 
presented with the exception of cases in which essentially identi- 
cal diagrams resulted from the three techniques. Certain features 
of the computed diagrams were discussed with respect to physical 
or mathematical. interpretation , but no conclusions were drawn. 

INTRODUCTION 

Man's increasing use of continental shelf regions of the 
world has stimulated the development of analytical models to aid 
in the study of the physical processes affecting offshore and 
coastal activities. One such analytical tool is the wave refrac- 
tion model, which describes the behavior of surface waves propa- 
gating through areas of irregular bottom topography. The Langley 



Research :Cehter (LaRC) and the Virginia Institute of Marine’ 
Science, (VIMS ) have, in a cooperative effort, developed a large- 
scale refraction model (ref. 1) for application. to the Virginian 
Sea, , a- region of' th.e mid-Atlantic continental • shelf from. Cape 
Hatteras .to Cape Heniopen. Computer resource requirements for 
applying the; model to a large geographical region. were signifi- 
cantly .reduced by the subsequent incorporation of a random-access 
technique for modular storage and retrieval of bathymetry ,d,ata 
(re.f . 2) . ' Recently , the model has been modified to . employ ,. oh 
optdon,^ three commonly used techniques fOr approximating Iqcal . 
ocean bottom topography. Two of these techniques , are- quadratic 
least squares . and cubic least squares , both of whiph pro-vide some 
smoothing of the discrete bathymetry data but do not insure con- 
tinuity in the depth or its derivatives. The third .technique , , 
constrained: bicubic interpolation , provides continuity without 
smoothing..^ Incorporation of all three techniques stemmed from a 
small-scale refraction study (ref. 3) in which it was concluded 
that;, the best ^approximation technique for use in future studies 
could riot .ibe determined without extensive comparison of computed 
results with; quantitative experimental data. In lieu of such a 
comparison,; refraction diagrams can be computed by using, each of 
the. ^bechriiques. Then more confidence can be placed in the:Cdm- 
puted Results ’ if no significant discrepancies exist among the 
three diagrams. , 

, - The purpose of this paper -is to present wave refraction 
diagrams computed by applying the LaRC-VIMS model to an area of 
the mid-Atlantic shielf called the Baltimore Canyon region after a 
prominent regional submarine geological feature. Extending from 
near . Wachapreague Inlet , Virginia, to near Manasquan Inlet, New 
Jersey, the region is of particular current interest with regard 
to offshore petroleum deposits and is one of a series of model 
regions which together provide broadrscale refraction modeling 
capability from Cape Hatteras to Montauk Point', Ne.w .York (ref. 4) 
An., immediate application of the computed diagrams is planned by 
the Philadelphia District and the Coastal Engineering' Research 



Center of the U.S. Army Corps of Engineers; these results wiil be 
used to furnish initial conditions for a fine-scale refraction 
study in the nearshore region off Ocean City, New Jersey. ■ • 

: DATA INPUT ' ’ 

The bathymetry data array for the Baltimore Canyon • regian-' was ' 
developed by first constructing a special transverse-Mercator map 
projection with a central meridian at 75® west, similar, to; the one 
constructed for the Virginian Sea region. (See ref ... l-pfor .;de tails .) 
The projection was then overlaid with a 0. 5-nauti.cal-^raile r'sjquare 
grid pattern extending for 150 nautical miles from north-' td south 
and for 127 nautical miles from west to east. Depth values ht the- 
76 200 nodes (300 rows by 254 columns) of the square grid. pattern 
were determined by using sounding charts and- other data* i*h?^a . man- -- 
ner analogous to that used in developing. ;the Virginian 'Sesi-:Tbathym- 
etry data grid (ref. 1). A three-dimensional computer-drawn .-plot "■ 
of the Baltimore Canyon region bathymetry data grid (with a verti- 
cal exaggeration of approximately 300 'to 1) is shown in figure 1.-' 
The regional array was then divided into a series of Smaller mod- '' 
ules for use with the random-access data storage and retrieval ' ' 
technique described in reference 2. The module size was selected"'’ 
as 42 rows by 32 columns to minimize computer cost,, resulting in 
a series of 72 overlapped bathymetry data modules. (See ref.- 2' 
for sequencing and overlapping details.) As a final adjustment to 
the input bathymetry data, the 'program input parameter for. tide 
was set so that the input depth values would correspond to mean 
low water . 

A range of wave direction and period was desired which would’ 
provide meaningful refraction diagrams while being realistic from- 
the standpoint of normally- occurring real-life wave conditions. 

With these requirements in mind, the selected input wave direc- 
tions ranged from 0° (from due north) through 90® (from due east) 
to 160® (from due south). Input directions from the west were not 
deemed significant for this study because of the limiting fetch for 
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physical wave generation. Within the selected direction range, 
input .directions were chosen at 22.5° increments (i.e., N, NNE, 

NE, etc. ) as was done in the study presented in reference 1.. 

Input wave periods were selected in 2-second increments from 6 sec- 
onds ta 16 seconds. On- the basis of results of ^ reference 1, waves 
of ' periods . less than 6 seconds were expected to undergo very little, 
if' any, refraction; thus, refraction-:diagrams-"for these periods 


would be of little meaning. It was also felt that waves of peri- 
ods greater than- 16 seconds would occur so infrequently in nature 
that their inclusion in this study was not justified. As a final 
input, the initial wave height was fixed arbitrarily at 0.3 meter 
for all cases since, according to linear wave theory, wave' height 
does not affect the refraction pattern for a given combination of 
wave period and direction. 

WAVE, REFRACTION DIAGRAMS ' 


Wave ref raction; diagrams • for' the Baltimore Canyon region were 
computed for the input range of wave period and direction • by using 
quadratic least squares, cubic least squares, and constrained 
bicubic interpolation techniques for approximating local bottom 
topography. The initial separation distance between waVe rays was 
selected as 3 nautical miles to provide adequate overall resolu- 
tion without cluttering the diagrams. Computations were required 
to begin in deep water, which, for the purposes of this paper, was 
considered 'to be a water depth greater than one-fourth the initial 
wavelength for the particular wave period being considered. (See 
ref. 1 for discussion.) Tick marks were drawn along each computed 
ray at points equally spaced in time, with the time increment 
varied with wave period to provide relatively uniform spatial reso- 
lution among- cases of different periods. By connecting correspond- 
ing tick marks o.n adjacent rays so that the connecting curve is 
orthogonal to each, ray, wave crest patterns can be derived. 

Computed- refraction diagrams are presented in figures 2 to 52 
in the sequence given in table I. Diagrams computed for a wave 
direction’ of 0° and periods of 12, 14, and 16 seconds have been 
omitted because only one or two rays met the deep-water condition 
for initiating computations. Unless noted otherwise, each figure;, 
consists of three parts: part (a) presents the diagram computed 

by using the quadratic least squares technique for approximating , 
bottom topography; part (b) presents the diagram computed by using 
the cubic least squares technique; and part (c) presents the dia- 
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gram computed by using the constrained bicubic interpolation tech- 
nique. When essentially identical diagrams resulted from the 
three techniques, only the one computed by using the quadratic 
least squares technique is presented. 

Some general comments are in order as an aid to potential 
users in interpreting the various diagrams. Since the tick marks 
along each ray denote equal time spacing, an increase in the spa- 
tial density of tick marks as a ray approaches the shore indicates 
decreasing wave phase speed. Convergence of a group of -neigh- 
boring rays indicates an area of relatively high wave energy,; 
whereas a sparsity of rays indicates relatively low energy. On 
occasion, individual rays stop abruptly in an area of relatively 
deep water (ray 23 in fig. 20(a), for example). Such ray: termi- 
nations indicate failure to converge on a proper value fo;r local 
ray curvature or calculation of a physically unrealistic value 
for wave height and, as such, should be considered of no physical 
significance. In several instances, individual rays undergo 
abrupt changes in direction relative to the mean direction of 
neighboring rays (ray 40 in fig. 32(c), for example). These 
abrupt changes are induced by sharp gradients in the local bottom 
topography, which are influenced, in turn, by the technique 
selected for approximating the topography within the finite-mesh 
bathymetry grid. Since use of the constrained bicubic interpola- 
tion technique results in a surface which passes directly through 
the bathymetry data points, sharp topography gradients tend to 
occur more often with this technique than with the two least 
squares techniques, both of which provide some smoothing of the 
data. 
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TABLE I.- SEQUENCE OF WAVE REFRACTION DIAGRAMS 


Figure 

Wave direction, 

Wave period, 


Figure 

Wave direction. 

Wave period. 


a, degrees 

T, seconds 



a , degrees 

T, seconds 

2 ; 

0 

a6 


28 

90 

16 

3 . 

0 

^8 


29 

112.5 

6 

4 

0 

10 


30 

112.5 

8 

5 . 

22.5 

®6 


. 31 

112.5 

10 

6 

22.5 

8 


32 

112.5 

12- ■ 

7 

■ . 22.5 

10 


33 

112.5 

14 

8 

22. 5 

12 


34 

112.5 

16 , 

9 

22,5 

14 


.• 35 

135 

6 

10 

22.5 

16 


36 

135 

8 

11 

45 

6 


37 

135 

12 

12 

45 

8 


38 

135 

13 

45 

10 


39 

135 

14 

14 

■ 45- 

12 


40 

135 

16 

15 

: 45 ' 

14 


41 

157.5 

6 

16 

45 

1 6 


42 

157.5 

8 

17 

'..67.5 

6 


43 

157.5 

10 

18 

,67.5 ,, 

8 


44 

157.5 

,12 

19 

,67.5 

10 


45 

157.5 

14 

20 

67.5 

12 


46 

157.5 

16 

21 

’ ■ ■ 67.5 

14 


47 

180 

6 

22 

67.5 

16 


48 

180 

8 

23 

90 

6 


49 

180 

10 

24 

90 

8 


50 

180, 

12 

25 

, ^ 90 

10 


51 

180 

14 ' 

26 

■ : 90' 

12 


52 

180 

16 

27 

' 90 

14 
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®Only the diagram computed by using the quadratic least squares 
technique is given because no resolvable difference existed among dia- 
grams computed by using the three approximation techniques. 
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(b) Bottom topography approximated by -cubic 
least squares technique. 

Figure 7.- Continued. . 
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X-axis, nautical miles 


Y-axis, nautical miles 
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(b)' Bottom topography approxima 
least squares techniq 
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Figure 9.- Continued 
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•(c) Bottom topography approximated 
;. bicubic interpolation tec 

Fieiure 9.- Concluded 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 10,-' Concluded, 
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X-axis, nautical miles 
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(a) Bottom topography approximated by quadratic •. 

•- - .. least squares technique . • . '• 

Figure 12 . ~ -Wave -refraction diagrams. • a = 45°..;'. ' *4 e d ends. 
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(a) Bottom topography approximated 
least squares technique 
: II;- Wave refraction diagrams. a = 4 





X-axis, nautical miles 







X-axis, nautical miles 
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X-axis, nautical miles 





X-axis, nautical miles 







X.-axis, nautjcal miles 


Y-axis, nautical miles 


DELAWARE 


APE;MAV 






/ C/" 




>.o‘ / 


I 






/M 



(b.) Bottom topography approximated by cubic 
least squares technique . 

Figure 15.- Continued. 
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(a) Bottom topography app 
least squares 
17. - Wave refraction diagra 
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( c )■ Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 1 8.- Concluded. 






(a) Bot-tom topography approximated by quadratic 
least squares technique. 

Figure 19. - Wave ' refraction diagrams, a = 67,5°; T = 10 seconds. 
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figure 19 .. \ 
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X-axis, nautical miles 
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(b) Bottom topography approximated by cubic 
least squares technique. 

Figure 21 -Continued . 





(c) Bottom topography approximated by constrained 
bicubic int'erpolatibn technique.' •“ .t-- 

• Figure; 21. -^Concluded. 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 23.- Wave refraction diagrams, a = 90°; T = 6 seconds 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 23 .- Concluded. 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 24.- Wave refraction diagrams. a = 90°; 
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T = 8 seconds. 







X-axis, nautical miles 



(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 24.- Concluded. 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 25.- Concluded. 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 26.- Wave refraction diagrams, a = 90°; T = 12 seconds. 
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X-axis, nautical miles 








X-axis, nautical miles 






(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 27.- Wave refraction diagrams, a = 9QOj T = 


seconds 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 27.- Concluded. 









(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 28.- Wave refraction diagrams, a = 90°; T = 16 seconds. 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 30.- Wave refraction diagrams. a = 112.5°; 
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T = 8 seconds. 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 31.- Wave refraction diagrams, o = 112.5°; T = 10 seconds 
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X-axis, nautical miles 








(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 32.- Wave refraction diagrams, a = 112.5°; T = 12 seconds. 
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X-axis, nautical miles 







(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 32,- Concluded. 






■(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure -Wave refraction diagrams. a = 112.5°; T = 14 seconds. 
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X-axis, nautical miles 






X-axis, nautical miles 







(a) Bottom topography approximated by quadratic 
least squares technique. 


Figure 3 ^ 1 .- Wave refraction diagrams, 
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T = 16 second 





X-axis, nautical miles 






X-axis, nautical, miles 













(b) Bottom topography approximat 
least squares techi'niqu 
Figure. 35.- Continued 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 36.- Wave refraction diagrams, a = 135°; T = 8 seconds. 
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(c) Bottom topography approximated by constrained- 
bicubic interpolation technique. 

Figure 36.- Concluded. 






(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 37.- Wave refraction diagrams. a = 135°; T = 10 seconds. 
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(b) Bottom topography approxiraat 
least squares techniqu 
Figure 37.- Continued 
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(a) Bottom topography approximated by quadratic 

least squares technique. - 

Figure 38.- Wave refraction diagrams, a = 135°; T = 12 seconds. 
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X-axis, nautical miles 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 39.- Wave refraction diagrams, a = 135°; T = 1^ seconds. 
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X-axis, nautical miles 













(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 40.- Wave refraction diagrams. a = 135°; T = 16 seconds. 
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(b) Bottom topography approximat 
least squares techniqu 
Figure 40.- Continued 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 40.- Concluded. 






(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 41.- Wave refraction diagrams, a = 157.5°; T = 6 seconds. 
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(b) Bottom topography approximat 
least squares techniqu 
Ficure 41.- Continued 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. ■ 

Figure 41.- Concluded. 







(b) Bottom topography approxiraat 
least squares techniqu 
Figure M2.- Continued 
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(a) Bottom topography approximated by quadratic 
least' squares technique. 

Figure 43.- Wave refraction diagrams, a = 157.5°;- T = 10 seconds 
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Figure 43.- Concluded 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

- Figure 44.- Concluded. 








X-axis, nautical miles 



; ,(a) Bottom topography approximated by quadratic 

' least squares technique. 

Figure; Wave refraction diagrams. a = 157.5°; T = 1 seconds 
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(a) Bottom topography approximated by quadratic 
least squares technique. 

Figure 46.- Wave refraction diagrams, a = 157.5°; T = 16' seconds 
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(c) Bottom topography approximated by constrained 
bicubic interpolation technique. 

Figure 46.- Concluded. 
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(b)' Bottom topography approximated by cubic 
least squares technique. 

Figure 48.- Continued. 
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(b) Bottom topography approximated by cubic 
least squares technique. 

Figure 51.- Continued. 
















(b) Bottom topography approximated by cubic 
least squares technique. 

Figure 52,- Continued. 
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edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTIUZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




